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^ , Abstract 

\_ j , By introducing generalized Backlund Transformations depending on arbi- 

Sr ■ trary functions, wave and localized soliton solutions of the Davey-Stewartson 

^"^ ' equations are generated. Moreover explicit soliton solutions of the Hamilto- 

J.—S , nian DSI and DSIII equations are obtained. 

1 Introduction 

O, 

C/3 _ Many efforts have been made in the last years to extend the soliton theory to the 

non Hnear evolution equations (NLEEs) in 2 + 1 (two spatial and one temporal) 
C^ I dimensions. The Spectral Transform (ST), which is the principal mathematical 

Oh' tool of the soliton theory, has been extended to dispersive NLEEs in 2 + 1 dimen- 

sions in the first '80 years, but it was generally admitted the lack of two dimen- 
. . sional localized solitons. Only recently, in 1988, it has been discovered by Boiti- 

rN I Leon-Martina- Pempinelli H] that all the equations in the hierarchy related to the 

jrt ' Zakharov-Shabat (ZS) hyperbolic spectral problem in the plane have exponentially 

localized soliton solutions. The most representative equation in the hierarchy is 
the Davey-Stewartson I (DSI) equation, which provides a two dimensional general- 
ization of the non linear Schrodinger (NLS) equation. These 2-diniensional soliton 
solutions (also called dromions ^j) display a richer phenomenology than in 1 + 
1 dimensions. The scattering of the solitons can be inelastic and they can change 
shape and also exchange mass. They can also simulate inelastic scattering processes 
of quantum particles as creation and annihilation, fusion and fission, and interac- 



tion with virtual particles 
method proposed by Sabatier 



Moreover, as a relevant application of the spectral 
5| and further developed by Boiti-PempincUi-Sabatier 
^, it has been stated that an additional NLEE, called Davey-Stewartson III (DSIII) 
equation admits also localized soliton solutions with properties similar to those of 
the DSI equation 0|. Localized soliton solutions for the DSI equation were first 
found by using gauge Backlund Transformations (BT) KM- Later these solutions 
were also obtained by means of the Inverse Spectral Transform (1ST) 13,0,0,0 and 
direct methods |iq,0,0|. Now we want to reobtain them by using a generalized 
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gauge BT pJ . The advantage of this method is that in such a way it is possible to 
obtain soliton solutions for DSI and DSIII equations in the so-called Hamiltonian 
case. For this very interesting case we have not presently at our disposal the ST or 
the dressing method, so to get explicit solutions we must generalize the Backlund 
gauge in order to include the special form of the boundaries in the Hamiltonian 
case. 

Finally let us note that a review on multidimensional localized solitons and more 
references are given in [ [l4| . 

2 Davey-Stewartson I and III equations 

Let us write the DSI equation in his more general two component form in charac- 
teristic coordinates 

iQt + (yz{Quu + Qvv) + [A,Q]^0 (2.1) 

where 

/ g(u,«,t) \ ( A^^){u,v,t) \ 

^'\r{u,v,t) )' \ A(2)(u,«,t) j ^^-^^ 

u and V are the characteristic coordinates u = x + y,v = x — y. 
The equation DSI is compatible with the reduction 

q = sr, e = ± (2.4) 

(where r denotes the complex conjugate of r), which furnishes the so-called re- 
duced DSI equation, describing physical situations as in hydrodynamics or in plasma 
physics. 

The equation DSI can be obtained m, m as the compatibility condition between 



two Lax operators Ti and T2 which commute in the "weak" sense |15 

Tiip^O, [Ti,T2]ip^0. (2.5) 

Ti is the ZS hyperbolic operator in the plane 

T,^ = (2 f ^'' M + 0) ^ - (2-6) 



and T2 has the form 

2T,^ ^ |^a, + d',-dl + A+(^ ° "J" ) I (^ = -k'^a,. (2.7) 

Let us remember that the boundary conditions of the auxiliary field A can be 
arbitrarily chosen. The equation DSI in his standard version, i.e. with the boundary 
written as 

A^'\u,v,t) = -- du'{Q'), + a^^\v,t) (2.8) 

■^J -00 

A^'-\u,v,t)^\[ dv'{Q'')^ + a!i^\u,t), (2.9) 



admits localized soliton solutions 



The one-soliton solution has the form 



Q = -j^>^iVe" 



r = --^A^/pe 



where 



D = 27(coshCi + cosh ^2) + e^' 

6 = -^J'IU - Xiv + 2(Ai{A/ + ^J,R^J.I)t 

6 = M/" - -^/w + 2(AijA/ - i-iRPi)t 

= ^RU + Xrv + {Xj - Xr + nj - HR)t 



1 



7 



rVP- 



(2.10) 



(2.11) 
(2.12) 
(2.13) 
(2.14) 

(2.15) 



The complex parameters A = A_r + iXi, /i — fj,R + ifii are the discrete eigenvalue of 
the associated Zakharov-Shabat spectral problem and p, rj are arbitrary complex 
constants satisfying the conditions 7 £ IR and 7(1 + 7) > 0. 
Another relevant choice for the boundaries can be done 



^ \^ — 00 ^ + 00/ 



(2.16) 



A(2) = 



+ 00 



(2), 



dv'{QX + K"{u,t) 



In this case one can introduce the Hamiltonian 



1 



H = dudv 
and the canonical Poisson bracket 



,(2) 



r{dl + dl)q - -qr{dud-' + d,d-')qr + (<^ - A};')qr 



{F,G} ^iffdudv 



5F SG SF SG 
Sq Sr Sr Sq 



where q and r are the conjugate variables. Then the equations of motion 

qt = {q,H}, n = {r,H} 



(2.17) 



(2.18) 



(2.19) 



yield the DSI equation. It has been shown that the DSI equation for Aq = Aq = 
is completely integrable in the hamiltonian sense 11^, ll^, so this case is called 
Hamiltonian. Quantum extensions have also been done. However the problem of 
defining a ST is completely open. 

There is another NLEE that can be associated to ZS hyperbolic spectral operator 
in the plane and admits localized soliton solutions, the so-called DSIII eq uation 
[bHI ■ Indeed the hamiltonian version of this equation appears already in [|l8| and a 
bihamiltonian version in [n9|. The DSIII equation has the form 




[A, Q] = 
(Q^), 



(2.20) 
(2.21) 



Also this equation is compatible with the reduction q — er. It can be obtained as 
the "weak" compatibility condition between two operator Ti and T2. Ti is given, 



as before, by (2.6). T2 in this case takes the form 



:^2 , C|2 , 1 , / Qu \ [ ,2, 



2T2^ ^ i^zdt + dt,+dt + A+l^-^ j j ^ = ~^^- ^^-^^^ 

Also the DSIII equation in his standard version, i.e. with the boundaries chosen as 

1 /"" 

A'-'Hu,v,t)^-- du'{Q\,+a^o'Hv,t) (2.23) 



2 



00 



A^^){u,v,t) = -\r dv'{QX + al^\u,t), (2.24) 

admits localized soliton solutions |^, , that have the same shape of those of the 
DSI equation, but different time evolution. The one-soliton solution has the form 

q = -^Xive'', r = --^A^/pe"^^ (2.25) 

where 

D = 27(coshCi + cosh^2) + e^^ (2.26) 

Ci = -/"/" - ^/« + 2(A7?A/ - ^J.R^^I)t (2.27) 

6 = fJ-iu - Xiv + 2{XrXi + iiBPi)t (2.28) 

e^^iRU + XRV + {x]~xl~^l]^^l]^)t (2.29) 

7 = \np- (2.30) 

As in the DSI case, the complex parameters A = A;? + iA/, /i = ^r + i/i/ are the 
discrete eigenvalue of the associated Zakharov-Shabat spectral problem and p, rj 
are arbitrary complex constants satisfying the conditions 7 S IR and 7(1 + 7) > 0. 
Moreover the choice for the boundaries 

Ad) = -\{r + /" ) du' {Q\ + Ai'\v, t) (2.31) 

furnishes the Hamiltonian case. If we introduce the Hamiltonian 



H = dudv 



r{-dl + dl)q + \qr{dud-^ - d,dz')qr + (A^ - A^^^r 



(2.32) 



the equations of motion qt = {9, H} and rt ~ {r, H} yield the DSIII equation. 

Finally let us make some com men ts. The DSI and DSIII equations belong to the 
same principal spectral problem (|2.6| ) , so the direct and the inverse problems are the 
same for these two equations. On the contrary, they have different auxiliary spectral 
problems, so the evolution of the ST and the dispersion relation are different. Also, 
as we have seen, there is a simple modification for the boundary conditions. The 
conclusion is that all the results known for the DSI equation can be easily extended 
to the DSIII equation . 



3 Solutions via Backlund Transformations 

Given a solution Q of the DSI (or DSIII) equation, we want to generate a new 
solution Q' of the same equation by introducing a convenient gauge operator B 
which transforms by means of 

^' = B{Q\Q)i: (3.1) 

the matrix solution ip of the principal spectral problem 

Ti(Q)^ = (3.2) 

for Q to the matrix solution ip' of the same spectral problem for Q' 

T^{Q'W = Q. (3.3) 

It is easy to verify that if B satisfies 

T^{Q')B{Q',Q)-B{Q',Q)T^[Q) = Q (3.4) 

T2[Q')B{Q', Q) - B{Q', Q)T2{Q) = 0, (3.5) 

then Ti{Q') and T2{Q') satisfy the same compatibility condition and therefore Q' 
satisfies the same equation as Q. The above equations furnish, respectively, the 
so-called space and time component of the BT. 

We are interested in the most general Backlund gauge as polynomial of the first 
order in dy 

B{Q',Q) = ady + BoiQ',Q) (3.6) 



with a a constant diagonal matrix and Bq a matrix. By inserting it in ( p.4[ ) , where 
Ti is the ZS hyperbolic spectral operator (p.6|), we get 

B{Q', Q) = ady - \<Jz{Q'a - aQ) - ^a^aliQ'^ ~ Q^) + P (3.7) 

and the space component of the Backlund transformation 

Q' [(3 - ^a<j3 1(Q'2 - Q2)] _ [^ _ laas I{Q'^ ~ Q^)]Q 

~\oz {Q'a - aQ), - i {Q'a + aQ)y = 0. (3.8) 

The matrix operator X is defined by 

I ^ [d^ + asdy)-' (3.9) 

and the diagonal matrix /3 is subjected to the constraint 

{d,+a3dy)P = 0, (3.10) 

i.e. it is of the form 

where pi and P2 are arbitrary functions. Note that for /3 a space dependence is 
allowed in contrast with the 1 + 1 dim case where /3 is the constant of integration. 
This freedom will be used for getting soliton solutions for Hamiltonian DSI and 



DSIII equations. By inserting (3.7) in (3.5) we get the time component of the 



BT, which is equivalent to the DSI (or DSIII) equation for Q' plus two additional 
equations that can be used for determining the auxiliary field A' and the admissible 



/3's. For details for the DSI equation see [13ll4|. The final result is that a and /3 
can be written as 

in the Backlund gauge B = B{Q,Q';X,fi). Note that A and /i generalize the pa- 
rameters of the BT in 1 + 1 dimensional case. 

For A and /z complex constants, the BT furnishes the localized soliton solutions 
for DSI (or DSIII) equation for the choice Q = 0, A = O[0,|. 

But to obtain soliton solutions for the Hamiltonian DSI (or DSIII) equation, we 
must use 

1. the space dependence of A and /i 

2. A can be ^ when Q = 0. 

In the case where A and /i have a space dependence, the calculations are very 
complicated and some constraint equations appear. We will give here only the 
principal results. Details for the DSI case can be found in [|l3|. We will restrict 
here at this equation. Results for the DSIII equation can be obtained with some 
modifications from those of the DSI equation. 

The starting solution is Q = and A = diag(AgQ^ (w, i), Aqq (u, t)) where A^^ are 
arbitrary boundary values. If we choose Aqq real and moving with constant speed 

4'oHu,t)=4o(" + 2^0, (3.13) 

a new solution Q, A can be explicitly written 

yl(i) = A« + 2^2 log(£P + pry/4), 
A(2)=4^„)-2a2log(fP + pr;/4) 
with boundary values 

4'^ = 4o + d^ log(f 2? + ot/4)1„=-oo + 32 log(£:i? + w/4)|«=+oo, 
4'^ = 4'o' - dl log(f P + p77/4)|.=_^ - dl \ogiSV + pvmv=+oo 

and 

d ^ (t)v + {(j)^ - (l)o)t + So, e = 9u + (6^ - eo)t + eo, (3.17) 

where (Jjq, 0q, Sq, eg are real constants. For the reduction r = eg, the wronskian W 
becomes constant 

W{r],V)^2a, W{p,£) ^ ~2ea, a £ M. (3.18) 

The real functions P, ?/ and £, p satisfy the stationary Schrodinger equations 

Vvv + (^00 ~ Mv = 



(3.16) 



and 

l(2) 

,(2)":"r":' (3-20) 



£-uu — (^00 + ^0)'? -- 0, 
Pun ~ (^S' + ^0)P = 0. 



The localized soliton solution is reobtained with a special choice of the functions 2?, 
77, £", p and with Aqq = Aqq = 0. We are interested in the Hamiltonian case 

4'^ = 4'^ = (3.21) 

in the reduced case. Let us introduce two functions 779 (w + 20t) and po{u + 2dt) 
defined as follows 

A^iJ = 2dl log ,70, ^S' = -25^ log Po. (3.22) 

We have to solve a complicated nonlinear system of coupled equations with con- 
straints for 2?, ?7, ?7o and f , p, po- If we add the additional constraints 

liin ^ = -2(pi±p2), (3.23) 

lim -^ = -2(771 ±772) (3.24) 

(pi and ?7i are real constants to be determined), the system decouples and we obtain 
the solution 

^ '2ap2V2PoVo exp[-7(e + S)] ^5) 

sinh a sinh (3 + (pi sinh a + a' p2 cosh a) (771 sinh /3 + ct"?72 cosh (3) 

where (cr')^ = 1; i'^")'^ = 1 ^^'^ '^^ P ^^^ related to T> and £ respectively. 770, a, po, 
/3 are determined by 

dim + a^p|?7§ - (/io»7o = 0, 9t,a = a'ap2r]o- 
dlpo + o^tiIpI - Oopo = 0, du(3 = -ao(j"aT]2pl 

with the consistency conditions 



(3.26) 



lim (pi +(T'p2C0tha) = -(?7i ±772) \ (3.27) 

lim (r/i+CT"772Coth/3) = -(/7i±/72)"^ (3.28) 

u — *±oo 

The above ordinary differential equations for r/o and po can be explicitly integrated 
in terms of elementary or classical trascendental functions and it is easy to verify 
the consistency conditions. 

Let us consider two cases in particular 



1. 



We get 



where 



d^r,o = 0, duPo = 0, 00 = A2 > 0, 0o^pI> 0. (3.29) 

(3.30) 



cosh^ 



£_ = po{u + 26t — ?io) - Ao(7; + 2(j)t — vq), 
e^9u+{0^ -pl)t + eo, 

S = (l)v + ((f>^-Xl)t + So, (3.31) 

so q in this case is a wave soliton solution. 



2. 

a„77o ^ 0, duPQ ^0, 00 < 0, 00 < 0. (3.32) 

In this case the solution can be exphcitly written p3[ in terms of the p, C, ct- 
Weierstrass ehiptic functions and describes an infinite wave with a periodically 
modulated amplitude. 
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